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Abstract: We observe that in SUSY models with non- universal GUT scale gaugino mass 
parameters, raising the GUT scale SU{2) gaugino mass IM2I from its unified value results 
in a smaller value of —m'jj^ at the weak scale. By the electroweak symmetry breaking con- 
ditions, this implies a reduced value of /x^ vis a vis models with gaugino mass unification. 
The lightest neutralino can then be mixed Higgsino dark matter with a relic density in 
agreement with the measured abundance of cold dark matter (DM). We explore the phe- 
nomenology of this high IM2I DM model. The spectrum is characterized by a very large 
wino mass and a concomitantly large splitting between left- and right- sfermion masses. In 
addition, the lighter chargino and three light neutralinos are relatively light with substan- 
tial higgsino components. The higgsino content of the LSP implies large rates for direct 
detection of neutralino dark matter, and enhanced rates for its indirect detection relative 
to mSUGRA. We find that experiments at the LHC should be able to discover SUSY over 
the portion of parameter space where rrig ~ 2350 — 2750 GeV, depending on the squark 
mass, while a 1 TeV electron-positron collider has a reach comparable to that of the LHC. 
The dilepton mass spectrum in multi-jet + £^i~ + E^^^^ events at the LHC will likely show 
more than one mass edge, while its shape should provide indirect evidence for the large 
higgsino content of the decaying neutralinos. 
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1. Introduction and framework 



There is overwhelming evidence showing that most of the matter in the Universe is not 
baryonic, but composed of a new massive stable (or at least very long-lived), weakly (or 
super-weakly) interacting particle not contained in the Standard Model (SM) of particle 
physics. Moreover, the mass density of this cold "dark matter" (DM) has been precisely 
measured: combining the results from the WMAP Collaboration with those from the Sloan 
Digital Sky Survey gives |0] 

njiuh^ = O.Ult^of^l {2a), (1.1) 

where = p/pc with pc the closure density of the Universe, and h is the scaled Hubble 
parameter, h = 0.73 it 0.04. While the mass density of DM is rather precisely known, the 
identity of DM remains a mystery. Like visible matter, the DM may well consist of several 



components, in which case the density in (1.1) yields an upper limit on the density of any 



single component. Supersymmetric (SUSY) models of particle physics with a conserved 
i?-parity include a stable, massive interacting particle, often the lightest neutralino Zi. 
Remarkably, the properties of the neutralino are just right to enable it to serve as thermally 
produced dark matter within the standard Big Bang cosmology if the SUSY mass scale is 
~ 100 GeV |, |.i 

A potential sticking point in the discussion of SUSY DM is that the non-observation 
of direct or indirect effects of SUSY are beginning to push MsusY beyond 100 GeV. As 
a result, the neutralino annihilation cross-section which is oc l/Mgygy is correspondingly 
reduced, and the calculated neutralino relic abundance is typically considerably larger ||5| 



than the value in Eq. (LI). For instance, in the mSUGRA model [^, almost all of parameter 
space is excluded by the precisely measured abundance, and only a few distinct regions 
where neutralino annihilation is enhanced survive: the nearly excluded bulk region with 
low masses already mentioned above Q, the stauQ or stop[^ co-annihilation regions, the 
hyperbolic branch/focus point (HB/FP) region at large mo||l^, where the Zi becomes 
mixed higgsino dark matter (with enhanced annihilation to W and Z particles via its 
higgsino component), or the A or h resonance annihilation (Higgs funnel) region [0, |l2|. 
Within the mSUGRA framework, each of these regions leads to characteristic patterns in 
collider signals at the soon-to-be-operational CERN Large Hadron Collider (LHC), and 
also to differences in signals in direct and indirect DM search experiments under way. 

We recognize, however, that these conclusions about the expected patterns are specific 
to the mSUGRA model, and can be obviated in simple extensions of this framework where 
new mechanisms may arise to match the predicted neutralino abundance to the measured 
CDM density. In previous works, in the scalar sector, non-universal soft masses for the 
different generations , or for the Higgs scalars have been considered |p^. In the gaugino 



^It is possible to invoke models where the DM results from late decays of heavy particles and so is not 
in thermal equilibrium, or to invoke non-standard Big Bang cosmologies |^] to obtain the observed value 
of the DM density. It is, nevertheless, appealing if we can account for the relic density data with minimal 
assumptions: i.e. via thermally produced relics within the standard Big Bang model. 
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sector, by abandoning gaugino mass universality at the GUT scale, a variety of new mecha- 
nisms emerge.^ As shown in Ref. |l^, 17], allowing the weak scale gaugino masses Mi ~ M2 



gives rise to mixed wino dark matter (MWDM), wherein ZiZi — > W'^W~ is enhanced in 
the early Universe. Alternatively, if GUT scale parameters are such that Mi ~ —M2 at the 
weak scale, then there is little bino-wino mixing, but neutralino annihilation is enhanced in 
the early Universe because of bino-wino co-annihilation p8|] (BWCA). Finally, if the SU (3) 
gaugino mass M3 <C Mi ~ M2 at the GUT scale, then the Higgs mass parameter m'jj^ is 
driven to less negative values at the weak scale, so that //^ ~ ~'^'h^ ^^^o small, resulting 
in mixed higgsino dark matter (MHDM) |ll9| , pO| ]. These models, wherein the composition 
of the neutralino is adjusted to get the correct dark matter abundance, are collectively 



dubbed "well-tempered neutralino" models] 21 1, with typical neutralino- nucleon scattering 
cross sections a{Zip) ^ 10~^ pb [22|, which is within an order of magnitude of the sensi- 
tivity of current experiments. It is instructive to note that each of these alternatives can 
arise in the top-down approach of string-inspired mixed moduli-anomaly mediated SUSY 



breaking (mirage unification) models] 23 1. It should also be noted that the low M3 frame- 



work with its concomitantly light squarks offers a novel possibility for getting agreement 



with (1.1) if the rate for ZiZi — > tt is unsuppressed because ti is relatively light [24|; the 



phenomenology of such a scenario has recently been detailed in Ref. |25| 



In this paper, we present a novel possibility, again based on non-universal gaugino mass 



parameters at the GUT scale, for obtaining agreement with (LI) via MHDM. We assume 
that the MSSM as the correct effective field theory valid between energy scales Q ~ M^eak 
and Q = Mqut ~ 2 x 10^^ GeV (for text book accounts, see Ref. |2|, |7|). We adopt the 
usual universality of scalar mass parameters and trilinear scalar couplings but take GUT 
scale boundary conditions in the gaugino sector of the form M2 ^ Mi ~ M^. We then 
inspect the evolution of the soft term rn'jj^ whose 1-loop RGE is given by 

where t = logQ^, ft is the top quark Yukawa coupling, Xt = fnq^ + m? -|- m'jj^ + and 
S = m'jj^ — fn^j^ + Tr (hiq — — 2m^ + + m|;) . Usually, the ffXt term overcomes the 
upward push from the gauge-gaugino terms (proportional to the gaugino mass parameters) 
and drives to lower values as Q is reduced, and ultimately to negative values, the 

celebrated radiative electroweak symmetry breaking (REWSB) mechanism]^]. In the case 
where M2 is very large at the GUT scale, the gaugino terms initially win resulting in an 
upwards push at the start of the fn^u^ evolution. The large M2 also increases the various 
left- scalar soft masses to initially large values, thus enhancing the magnitude of Xt which 
results in an increased downward push of the ffXt term. The resulting value of rn?^^ (weak) 
depends on the value of M2(GUT); by adjusting the latter we can arrange things so that 
by the time the weak scale is reached, the (incomplete) cancellation between the upwards 
and downwards push results in a negative value of rn\j^ that is significantly smaller in 



^Non-universal gaugino masses can be accommodated in SUSY GUT if the auxiliary field that breaks 
supersymmetry also breaks the GUT symmetry ]l5(]. 
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magnitude than in models where GUT scale gaugino mass parameters have a common 
value. The weak scale value of fi^ (at tree-level) is then obtained from the weak scale 
parameters of the Higgs sector via the EWSB relation, 

^ " (tan2/3-l) ^^-"^^ 

We see that if ^ and for moderate to large values of tan/3, ~ —m'jj^. 

Thus, the small value results in a smaller parameter, and a correspondingly 

larger higgsino component of the lightest neutralino Zi. As discussed above, MHDM can 
easily be compatible with the observed relic density because the neutralino annihilation 
into WW, ZZ and Zh pairs is enhanced on account of the higgsino content of Zi. Since 
the large value of the wino mass underlies the root of this scenario, we will call it the high 
M2 dark matter (HM2DM) model. 

The situation is illustrated in Fig. ||. The model is specified by the mSUGRA parameter 
set augmented by one additional parameter, M2: 

mo, 5 ^2, Aq, tan/3 and sign(^) , (1-4) 

where Mi = M3 = mi/2, but M2 is allowed to be free (with either sign). We take nit = 
171.4 GeV, in accord with recent mass determinations (2^. We illustrate in Fig. |l]a) the 

2 \ /l^2 I o;„v,/^2 \ /|^2 



evolution of sign(mj^^)y |m|^^| and sign{mjj^)^\mjjj as a function of scale Q in the 
mSUGRA model (solid) for mo = 300 GeV, mi/2 = 300 GeV, ^0 = 0, tan /? = 10, /x > 0. 



m„ =300aeV, m,,^ =300GeV, (anp=ro, A„ =0, ^ >0, m, =171.4GeV 
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Figure 1: Evolution of a) the soft SUSY breaking Higgs mass parameters sign(m|^^)y^|m|f | 

and sign(m|^^ ) ^ \ mf^^ \ , and b) the gaugino mass parameters, as a function of the energy scale Q 
in the mSUGRA model (solid) for toq = 300 GeV, toi/2 = 300 GeV, Aq = 0, tan/3 = 10, ^ > 
and mt = 171.4 GeV, and for the HM2DM model (dashes). The model parameters adopted for 
HM2DM are the same as in the mSUGRA case except that M2 = 3mi/2 at Q = Mqut- 



a (GeV) 
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The same running mass parameters are shown for HM2DM for the same parameters as in 
the mSUGRA case except that we now take M2 = 3mi/2 at Q = Mqut (dashes). In the 
mSUGRA case, fnj^^ evolves from a positive GUT scale value to a large negative value at 
Q = M„eak) resulting in a large |^| parameter. In the case of HM2DM, however, the large 
value of M2 initially causes rn^u^ to evolve upwards, but ultimately, the ffXt term wins 
out and rn?fj^ evolves to a not-as-large negative value and electroweak symmetry is broken. 
The smaller value of —m'jj^ at the weak scale leads, of course, to a correspondingly smaller 
value of 1^1 compared to the mSUGRA case. In Fig. we show the evolution of gaugino 
masses in mSUGRA and in the IIM2DM case. For mSUGRA, we are left with the gaugino 
masses at the weak scale in the well-known ratio of Mi : M2 ■ ~ 1 : 2 : 7 which implies 
(when is large) that m^^ : ^n-pf;^ : rrig ~ 1 : 2 : 7. In contrast, in the IIM2DM scenario, 
we have at the weak scale Mi <^ M3 ~ M2, so that we expect in general the SU{2) winos 
to be similar in mass to gluinos, and hence almost decoupled at the LHC.^ 




Figure 2: a) Evolution of the soft SUSY breaking squark mass parameters TOj^ and m^^ as a 
fmiction of scale Q in the mSUGRA model (solid) for toq = 300 GeV, mi/2 = 300 GeV, ^0 = 0, 
tan/3 = 10, /X > and mt = 171 A GeV. The same running mass parameters are shown for HM2DM 
model for the same parameters as in the mSUGRA case except that we take M2 = 3mi/2 at Mqut 
(dashes), h) Evolution of third generation soft SUSY breaking slepton mass parameters Wfj^ and 
TOf^ as a function of scale Q for the same case as in frame a). 

The large value of M2 in the IIM2DM model also influences the evolution of matter 
scalar mass parameters. In Fig. |2|a), we show the evolution of third generation squark 
mass soft parameters versus energy scale for the mSUGRA and the HM2DM models. In 
the mSUGRA case, we see the hierarchy m? > > generated by the large top 



^The reader may well wonder whether it is possible to obtain MHDM by increasing Mi instead of M2. 
This does not, however, appear to be possible because the bino mass required for this is so large that, though 
/i is indeed reduced, the lightest neutralino dominantly becomes a wino-higgsino mixture and annihilates 
too rapidly to saturate (1.1). 
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quark Yukawa coupling. However, in the HM2DM case, the large value of M2 causes left- 
matter scalars (which have gauge couplings to the heavy winos) to evolve to larger values, 
so we get large weak scale value of m^^ (= w-f,^), and also for first and second generation left- 
sfermions. In addition, the large m^-^ enters via the ffXt term in the RGE for m? , causing 
it to reduce at the weak scale. Thus, in HM2DM, we expect large L — R splitting in the 
squark sector, with mi^{= mj^^) > rn^^ > ^Ir- ^ result, in contrast to the mSUGRA 
framework, the lighter sbottom quark 61 is dominantly bn in the HM2DM scenario. 

In Fig. |2|6), we show the evolution of third generation slepton masses. In the mSUGRA 
case, there is only a small intra-generation splitting at the weak scale. In the HM2DM 
scenario, mf^ evolves as in mSUGRA, but m^^ gains a big enhancement from the large 
value of M2. Thus again, in the slepton sector, we expect large left-right splitting of slepton 
mass parameters in HM2DM, with left sleptons much heavier than the right sleptons. A 
similar behaviour is expected for the first two generations of sleptons. 

The upshot is that if a large M2 parameter is the underlying reason for MHDM in our 
Universe, then characteristic sparticle mass spectra and sparticle mixing patterns should 
emerge as a result. Our goal in this paper is to lay out the phenomenology of this framework. 
In Sec. ^ we explore the sparticle mass and mixing patterns, and delineate the allowed 
parameter space in the HM2DM scenario. In Sec. |3|, we examine the low energy constraints 
from 6 — > 37 and the SUSY contributions to the anomalous magnetic moment of the muon. 
In Sec. ^, we examine prospects for direct and indirect dark matter detection in the IIM2DM 
framework. In Sec. ^, we examine collider implications of the HM2DM scenario. We end 
in Sec. |6| with a summary of our results. 



2. Parameter space, relic density and mass spectrum 

As discussed in the last section, the HM2DM model is completely specified by the parameter 
set (Q), 

"^Oi '"T'i/2) ^h, ^0) tan/? and sign(/i), 

where we assume that Mi = M3 = mi 12 > at Q = Mqut, and where M2 can assume 
either sign. The assumed equality of Mi and M3 can be relaxed somewhat and our conclu- 
sions suffer little qualitative change so long as M2 ^ Mi. To calculate the sparticle mass 
spectrum, we use Isajet 7.76[pO|], which allows for the input of non-universal scalar and 
gaugino masses in gravity mediated SUSY breaking models where electroweak symmetry 



is broken radiatively. The relic density is evaluated via the IsaReD program] 31], which is 
part of the Isatools package. IsaReD evaluates all 2 — >■ 2 tree level neutralino annihilation 
and co-annihilation processes and implements relativistic thermal averaging in the relic 
density calculation. 

In the upper frame of Fig. ^, we show the neutralino relic density for the parameter 
space point mo = m.1/2 = 300 GeV, ^0 = 0; tan/3 = 10 and ^ > 0, versus the ratio 
r2 = M2(GUT)/mi/2- For r2 = 1 corresponding to the mSUGRA model, ^^i^'^ ~ 
that the point is strongly excluded because it yields too much dark matter. For r2 ~ 0.6, 
we arrive at the MWDM case explored in Ref. |l^], while for M2 ~ —0.5, we have BWCA, 
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explored in Ref. |18]. If we instead increase the magnitude of M2, then we find at r2 ~ 3 a 
match to the measured rehc abundance, with 0.7^ ^ 0.1. This is the case of HM2DM. As 
we increase r2 further, the rehc density starts to go back up. This is because the neutrahno 
mass falls below Mz and ultimately My/, so that (except neutralinos with high thermal 
energy) the processes ZiZi ZZ, WW become disallowed, and the total annihilation 
cross section is correspondingly reduced. Ultimately, the relic density again starts to drop 
because of the annihilation via off-shell Z. We note here that since the gaugino masses 
enter the soft SUSY breaking Higgs masses as their square, there is also a point with good 
relic density at r2 ~ —2.5. Again, the shoulder in the relic density curve marks where the 
annihilation to vector boson pairs becomes disallowed. 



trig =300GeV, /t?,/^ =300GeV, tan^=10, Ag =0, \i >0, =171.4GeV 




= M^/ 



Figure 3: The neutralino relic density ^cDAih^ (upper frame) and the higgsino component R^j 
of the lightest neutralino (lower frame) as a function of r2, for mo = 300 GeV, mi/2 = 300 GeV, 
^0 = 0, tan/? = 10, ^ > and rrit = 171.4 GeV. In the grey regions, the chargino mass falls below 
its lower limit of 103.5 GeV from LEP2 experiments, while the red hatched regions arc excluded 
either because the Z width becomes too large, or because electroweak symmetry is not correctly 



broken. The green region is where the relic density falls in the range (1.1) 



('l")2 ('1)2 

We also show the higgsino content of Zi, defined hy Rfj = v\ + V2 (in the 
notation of Ref. [^]) in the lower frame of Fig. |^. We see that over the bulk of the range 
of r2, the higgsino composition of Zi is quite low, since Zi is dominantly bino-like, or for 
< r2 < 0.6, a mixture of wino and bino. In the case of HM2DM with r2 — 3 or —2.5, Rp^ 
has risen to ~ 0.6, indicating mixed higgsino-bino dark matter. 
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In Fig. ^, we show the thermally averaged neutralino annihilation cross section times 
relative velocity, integrated from temperature x = T/m^^ = to freeze-out x = xp, versus 
r2 for the same parameter choices as in Fig. |3[ The inverse of this quantity enters the 
relic density calculation, so that a large integrated annihilation rate leads to a small relic 
density. For clarity, we display only positive values of r2. We see that while neutralino 
annihilation to lepton pairs via slepton exchange is dominant in the case of mSUGRA, 
when we move to the case of HM2DM, where the Zi is a mixed bino-higgsino state, then 
annihilation to WW, ZZ and Zh dominates, as is typical for mixed higgsino dark matter. 
As just discussed, these cross sections drop-off near the upper end of the range of r2 once 

falls below Mz or Mw. In this range, annihilation via s-channel Z (which has large 
couplings to the Zi pair on account of the large higgsino content of Zi) dominates, and 
ultimately becomes resonant so that the relic density drops below its observed value. 



nig =300GeV, m,^^ =300GeV, tan^=10, Ag =0, n >0, m, =171.4GeV 




Figure 4: Thermally averaged neutralino annihilation cross sections (calculated in the HM2DM 
model) times neutralino relative velocity, integrated from x ^ to xp versus r2, for the same 
parameters as in Fig. |^. 

In Fig. ^, we show the sparticle mass spectrum versus r2 for the same parameter 
choice as in Fig. |3|. At r2 = 1, we see a large mass gap — m^^ 100 GeV in the 
case of the mSUGRA model. As AI2 increases, the fi parameter decreases, and falls rapidly 
beyond r2 ~ 2. In the region of r2 ~ 2.5 — 3, the curves for m^^, m^^ and tti:^^, and for 
very large r2, in place of track the value, indicating that these particles are 
dominantly higgsino- like. We also see that mg^, and m- j all increase with increasing r2, 
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trig =300GeV, m^,^ =300GeV, tan^=10, Ag =0, n >0, m, =171.4GeV 




Figure 5: Various sparticle and Higgs boson masses and the /Lt parameter in the HM2DM model 
versus r2 for mo ~ 300 GeV, toi/2 = 300 GeV, An = 0, tan/3 = 10 and ^ > 0. 

owing to the upward push in their respective RGEs. The value of rriA, given at tree-level by 
~ m'jj^—m'j^^, also increases with r2 as can be seen from Fig. ||a. The right- squark and 
slepton masses, on the other hand, are roughly independent of r2, so that in HM2DM we 
expect a larger mass gap between L and R squarks and sleptons relative to the mSUGRA 
model. This also leads to the level-crossing in the 6-squark system that we mentioned 
earlier: for small values of r2 (including in the mSUGRA model) bi is dominantly bi, while 
for r2 ^ 2, bi becomes mostly The value of m^^ actually decreases with increasing M2, 
which is due in part to the diminishing value of m? as shown in Fig. 0a), and in part 
due to an increasingly negative weak scale value of At. For the most part, the figure is 
nearly symmetric between positive and negative values of M2 since, as mentioned above, 
the scalar SSB RGEs contain Ml, and not M2. The ^i,6,T parameter RGEs all have M2 
entering linearly, so that A term evolution is not symmetric between positive and negative 
fi. This gives rise to the unsymmetrical behavior, most noticeable in the m^^ curve. Of 
course, which becomes essentially IM2I in the HM2DM model- is also asymmetric 

since the value of \r2\ required to saturate the relic density is itself asymmetric between 
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positive and negative masses. This asymmetry could have an impact upon the accessibihty 
of W2 and Z4 at future electron-positron colliders. 

We have already discussed how the intra-generation mixing patterns are affected by 
the large value of IM2I: since a large value of results when IM2I is large, the lighter 

JL 

sfermions are dominantly fji within the HM2DM framework. This is confirmed in Fig. |6|, 
where we show the dependence of the sfermion mixing angle Of defined in Ref. |26| on r2 
for the same parameters as in Fig. ^. We see that the bi and fi are essentially b^ and f/j 
in the HM2DM model. 



irig =300GeV, m^^ =300GeV, tan^=10, Ag =0, \i >0, m, =171.4GeV 
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Figure 6: The sfermion mixing angle 9f as defined in Ref. [g6j for f = t, b and t versus Ti for 
mo = 300 GeV, mi/2 = 300 GcV, Aq = 0, tan/3 = 10 and ^ > 0. We see that for the HM2DM 
model, the 9f are all large so that the lighter states are all dominantly right sfermions. The various 
shadings are the same as in Fig. ||. 

In Table |l|, we list various sparticle masses and fi, along with expectations for the 
relic density, a^^^ , BF{b sj) and spin-independent direct DM detection cross section 
asi{Zip) for several different cases with all parameters other than M2(GUT) set as in 
Fig. 1^. The first case, mSUGRA, provides a benchmark for comparison with the HM2DM 
cases. For case HM2DM1, we take M2 = 900 GeV, to obtain ^^Ji^ ~ 0.1 as required 
by observation. As expected the gluino mass hardly changes in going from mSUGRA to 
HM2DM1, while the light charginos and neutralinos have all gotten much lighter in accord 
with the decreasing parameter. The light chargino also changes its character from being 
dominantly wino-like to dominantly higgsino-like. In fact, the mass gap m^^ — m^^, which 
was of order 100 GeV in mSUGRA, has dropped to ~ 50 GeV in IIM2DM1. This means 
the spoiler decay modes Z2 — > Zih and Z2 Z\Z will be closed for HM2DM, and a 
dilepton mass edge should be visible in collider events where Z2 is produced at large rates 
either directly or via gluino and squark cascade decays. The neutralino Zj, also has a large 
higgsino component and cannot be split very much from Z2; its leptonic decays, therefore, 
should also lead to a distinct mass edge in the dilepton mass spectrum. We also note 
that the left- squarks and sleptons have become 200-300 GeV heavier than in mSUGRA, 
while the masses of right- squarks and sleptons, as expected in the IIM2DM scenario, are 
essentially unchanged. The top squark t\ has become significantly lighter in the HM2DM1 
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parameter 


mSUGRA 


HM2DM1 


HM2DM2 


HM2DM3 


mo 


300 




300 




300 


300 


Ml 


300 




300 




300 


300 


M2 


300 




900 




-700 


-695 


M3 


300 




300 




300 


300 




385.1 




134.8 




136.5 


-144.4 




729.7 




736.4 




749.7 


749.5 




7z0.8 




901.8 




840.7 


838.6 




OZO.41: 




oy4.o 




533.0 


534.7 




DOD.O 








701.2 


700.7 








uuy.o 




559.9 


557.0 


men 


322.3 








QOI A 
oZi.4 


oZi.4 




411.7 




719.7 




0/0.1 


0(0.1 








1 ^fi 




1 QQ O 

ioo.y 


^AAA 
144.4 




41 9 




79^^ 1 
1 zo. 1 




cr 

583.2 


f on ^ 

580.7 




oyi.o 




1 fin 




170.1 


1 no A 

168.4 




zZU.D 




1 /1 Q 




136.2 


141.9 




1 1 n 

iiy.2 




94.8 




99.9 


108.6 


mA 


D/U.o 




(U. / 




565.0 


563.3 




^9Q 8 




fi7Q 8 




574.3 


572.7 


mh 


110.1 




111.9 




1 n7 K 


lU ( .0 




1.1 




0.10 




0.11 


0.12 


BF{b S7) 


3.0 X 10" 


-4 


2.3 X 10" 


-4 


3.3 X 10"^ 


4.0 X 10"^ 




12.1 X 10" 


-10 


3.1 X 10" 


10 


-7.4 X 10-1° 


7.0 X 10^1° 


crs/(Zip) 


2.1 X IQ-^ 


pb 


3.4 X 10"*^ 


pb 


2.5 X 10"^ pb 


3.2 X 10"^ pb 




0.05 




0.40 




0.37 


0.32 



Table 1: Input parameters and resultant sparticle masses in GeV units together with the predicted 
neutrahno rehc density, direct LSP detection scattering cross section from a proton, B{b 57) and 
Aa^, the SUSY contribution to the anomalous magnetic moment of the muon, for mSUGRA and 
three HM2DM scenarios. In each case, we fix ^0 = 0, tan/3 = 10 and nit — 171.4 GeV. 



case, which leads to a deviation of the branching fraction B{b ^ S'j) from its SM value. In 
contrast, the SUSY contribution to Aa^ is diminished, owing to the increased left smuon 
and sneutrino masses. The DM direct detection rate has increased to the ~ 10~^ pb 



level expected in models with MHDM [g2|]. In case HM2DM2, we take M2 = -700 GeV, 
which also gives the correct relic abundance. In this case, the ti is heavier than in case 
HM2DM1, so that BF{b — > ,37) is more closely in agreement with its measured value: 
BF{b 57) = (3.55 ± 0.26) x 10"^ from a combination of CLEO, Belle and BABAR 
data fs^ . However, the value of Aa?y^^ , which is proportional to "^'^■ft*^^ ^^'^ ^ , has turned 

SUSY 

negative, in contrast to the measured deviation, which is positive. This can be rectified 
by choosing in addition to M2 < 0, < 0, as in the HM2DM3 case in the last column 
of the table. This case now gives a positive contribution to Aa^, but also a deviation in 
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BF{b — > 57) which is now somewhat larger than the measured value. This latter case with 
opposite signs of ^ and Mi also gives a significantly lower direct DM detection cross section, 
due to negative interference between h- and -ff-mediated scattering amplitudes 

While our discussion up to now has been confined to particular values of mo, ' ' 

it should be clear that the method of raising IM2I to obtain MHDM in agreement with 

is quite general, although a different value of \r2\ will be needed for each point in 
rriQ vs. mi/2 parameter space. We have scanned over the mg vs. mi/2 parameter space for 
ylo = and tan/3 = 10 to extract the particular r2 value needed at each point to obtain 
the WMAP measured CDM density. The results are shown as contours of r2 in Fig. ^ for 
a) M2 > with /i > 0, and b) M2 < with n < 0. The red-shaded regions are theoretically 
excluded, due to lack of EWSB, due to a stau or (if r2 is very small) chargino LSP, or 
because the Z-width constraint is violated. The blue shaded regions are excluded because 

< 103.5 GeV, in contradiction with sparticle search limits from LEP2. The green 
shaded regions give /i^ < 0.13 in the mSUGRA model, so in these regions there is no 
need to dial M2 to large values. In frame a), the contours in r2 range from r2 ~ 3 in the 
low mo, mi/2 region, to r2 ~ 2 when nearing the HB/FP region (which already has \n\ 
suppression due to a large mo value). We note the appearance of a white region to the left 
of the mSUGRA stau co-annihilation region that is allowed in the HM2DM framework; 
in this case, since |^| is reduced, the value of m^^ ~ |^| and falls below m^^. There is 
also a region just below mi/2 ~ 0.3 TeV which turns out to be LEP2 excluded, while even 
smaller values of mi/2 ~ 0.2 TeV re-emerge as LEP2 allowed. The LEP2 excluded region 
at ini/2 ~ 225 — 290 GeV occurs because for mi/2 < 290 GeV, m^^ drops first below Mz, 
then below M^. This shuts off the ZiZi — > ZZ, WW annihilation modes, so that even 
larger M2 values are needed to drive |^| to even smaller values. The lower values then 
drive m^^ below the LEP2 search limit that requires m^^ > 103.5 GeV. For even lower 
mi/2 values, the allowed region opens up again because, for the smaller value of mi/2, the 
higgsino components of Zi allow for efficient annihilation through the off-shell s-channel 
Z exchange so need not be as small, and the chargino mass can be above the LEP2 
bound in the range of mi/2 ~150-225 GeV. 

The corresponding contours for negative M2 are shown in Fig. |^6), where we also take 
/i < to realize a positive value of Aa^^^^; the frame mainly differs from the positive M2 
case in that the range of \r2\ is somewhat lower than in frame a), presumably because of 
differences in m^^ and in the coupling of Zi to W^s and Z's. The LEP2-allowed region at 
mi/2 ~ 0.2 TeV from frame a) now also disappears (despite the smaller value of r2, fj. is 
now slightly smaller than in frame a, making it more difficult to evade the LEP 2 bound), 
except for the thin sliver where 2m^^ ~ mh where, because of the resonance enhancement, 
a raised value of M2 is not needed to saturate the measured relic density. 

Up to now, we have confined our discussion to a fixed value of tan/3 = 10. Our 
considerations also apply for other values of tan /3. Specifically, we have checked that for 
mo = mi/2 = 300 GeV, ^0 = and fj. > 0, that the relic density measurement is saturated 
for r2 : 2.8 — 3.4 for tan/? ~ 40; for yet larger values of tan/3, the required value of r2 
drops rapidly, and consistency with the upper bound on the relic density is possible in the 




Figure 7: Contours of r2 in the toq vs. toi/2 plane with tan/3 = 10, = for a) M2 > with 
/i > 0, and b) M2 < with /i < 0. Each point in these planes has r2 dialed to a high value such a 
value that f^^^/i^ ~ 0.11. The red and blue regions are excluded for reasons explained in the text. 



mSUGRA model once tan /3 ~ 46 because neutralino annihilation via s-channel A becomes 
large, while for the largest values of tan f3 co-annihilation via staus becomes dominant. 



3. b — > S7 and (g — 2] 



Now that we have established that any point in mo vs. space can be made dark- 

matter consistent by increasing IM2I, we delineate regions of parameter space where the 
recent measurements of the branching fraction BF{b — > 57) |^2| or of {g — 2)^ [35, 36| are 
consistent with predictions of the HM2DM model. 

3.1 BF{b S7) 

The branching fraction BF{b S7) is extremely interesting largely because amplitudes 
for supersymmetric contributions mediated by Witj and bH~^ loops are expected to be 



of similar size as the leading SM amplitude mediated by a tW loop |37|. The measured 
branching fraction, from a combination of CLEO, Belle and BABAR experiments p2| , 
is BF{b 57) = (3.55 ± 0.26) x 10"^, while the latest SM calculations find|3| BF{b 
57) = (3.29ib0.33) X 10~^. In view of the good agreement between the SM and experiment, 
any SUSY contribution to BF{b — > 57) should be somewhat suppressed, unless there are 
cancellations between different SUSY loops, or the summed SUSY contribution fortuitously 
turns out to be twice the SM amplitude but with the opposite sign. 

We evaluate BF{b 57) using the IsaBSG code [^], a part of the Isatools package. 
The dependence of the branching fraction on r2 is illustrated in Fig. ^ for our canonical 
point 1 from Table ||. We see that in the mSUGRA case, BF{b 57) is not far below 
its measured value. Dialing r2 to high positive values decreases the branching fraction 
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(recall that ii becomes lighter) and the discrepancy with experiment grows until we hit the 
green DM allowed region, where we find BF{b —>■ 57) ~ 2.3 x 10~^. For negative values 
of r2, BF[h sj) varies much less, and in the WMAP-favoured range is actually in close 
accord with its measured value. For the favoured negative sign of n, the branching fraction 
becomes a bit too large when r2 < 0. 



iVg =300GeV, /r7,/2 =300GeV, tan^=10, Ag =0, 11 >0, =171.4GeV 
4 




= M^/ rriy^ 



Figure 8: Branching fraction BF(b sj) versus M2 for mo — mi/2 
tan/3 = 10, and /i > 0. The various shadings arc the same as in Fig. ||. 



300 GeV, Ao = 0, 



While this particular point at high M2 > may seem somewhat discouraging, in Fig. ^ 
we show by the black contours the branching fraction BF(b 57) in the rriQ vs. mi/2 plane 
with ^0 = Oj tan /3 = 10 and a) M2 > with fi > 0, and b) M2 < with fi < 0, where at 
each point we have dialed IM2I to high values so that that ^^^h'^ saturates the measured 
value in ( |1.1| ). In frame a), we see that BF{b — > 57) is low only in the very low mo and 
mi/2 corner, but is not far from its measured value (considering theoretical uncertainties) 
for mo, mi/2 ~ 0.5 TeV. In contrast, in frame b), the measured value of BF{b sj) 
clearly disfavours small values of mo and mi/2) requiring these to be ~ 700 GeV. 

3.2 {g - 2)^ 

Current measurements of the muon anomalous magnetic moment show an apparent devia- 
tion from SM predictions. Combining QED, electroweak, hadronic (using e^e~ — > hadrons 
to evaluate hadronic loop contributions) and light-by-light contributions, and comparing 
against measurements from E821 at BNL, a positive deviation in 



(9-2)m 



of 



Aa,, 



exp 



22(10) X 10" 



-10 



(3.1) 



is reported in the Particle Data Book|35], i.e. a 2.2a effect."^ 

*More recent analyses|^^ report a larger discrepancy if only electron-positron data are used for the 
evaluation of the hadronic vacuum polarization contribution; the significance of the discrepancy is, however, 
reduced if tau decay data are used for this purpose. 
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Figure 9: Contours of branching fraction BF{b 57) x 10"' (black) and Aa^ x 10^" (purple) in 
the Too vs. mi/2 plane for Aq = 0, tan (3 — 10, and a) r2 > with /i > 0, and 6) r2 < with /i < 
where M2 has been dialed at every point to large values such that ^^i^^ ~ 0.1. 



One-loop diagrams with Wi — and Zj — in the loop would give super symmetric 
contributions to a^, perhaps accounting for the (rather weak, yet persistent) discrepancy 
with the SM. For our canonical point with rriQ = mx/2 = 300 GeV, we have checked 
that even though IM2I is large, the total chargino contribution to Aa^ dominates the total 
neutralino contribution even in the HM2DM model exactly as in the mSUGRA case. This 
is the case for both signs of M2. For M2 > the total neutralino contribution, though much 
smaller than the corresponding chargino contribution, is relatively larger in the HM2DM 
case as compared with mSUGRA. 

The purple curves in Fig. ^ are contours of Aa^^^^ x 10^^ in the mo vs. mi/2 plane. 
For M2 > shown in frame a), in the portion of the plane not strongly excluded by the 
6 ^ S7 constraint, we see that model, Aa^^^'^ is very small. The situation for M2 < is 
very similar as can be seen from frame h). While it appears that the HM2DM model will 
be strongly disfavoured if the muon magnetic moment discrepancy continues to persist, we 
should remember that (1) Aa^^^^ and BF{b — > S7) are both sensitive to tan /3, and (2) the 
latter is very sensitive to small flavour violations in the soft-SUSY breaking parameters 
which will not have any significant effect on direct searches for supersymmetry. 



4. Direct and indirect detection of neutralino CDM 

In this section, we explore the prospects for direct and indirect detection of neutralino 



dark matter within the HM2DM framework] 39]. We adopt the IsaReS code|2^ (a part of 
the Isatools package) for the computation of the direct detection rates and the DarkSUSY 
code [^], interfaced to Isajet, for the computation of the various indirect detection rates. 
Indirect detection rates are sensitive to the DM distribution in our galactic halo, larger 
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rates being obtained for a clumpy or cuspy halo distribution as compared with a smooth 
or less peaked distribution of the DM. We show our results for two halo profiles. The 
Adiabatically Contracted N03 Halo model[^], where the deepening of the gravitational 
potential wells caused by baryon in-fall leads to a higher concentration of DM in the center 
of the the Milky Way, gives higher detection rates, especially for gamma ray and anti- 
particle detection than smoother halo profiles. For comparison, we also show projections 
using the Burkert profile 1 43 1 where the central cusp in the DM halo is smoothed out by 



significant heating of cold particles.^ In our analysis, we consider signals from the following 
processes. 

1. Relic neutralinos in our galactic halo can scatter from nuclei in the material of un- 
derground cryogenic or noble liquid detectors designed to detect the resulting nu- 
clear recoil, leading to direct detection of the neutralino [^]. Although there is no 
positive signal to date, the most stringent upper limit on the scattering cross sec- 



tion comes from the XENON-10 collaboration |47|, which obtained an upper limit 
a{Zip) ~ 8 X 10~^ pb for m^^ ~ 100 GeV, corresponding to the expected neutralino 
mass in the HM2DM model for our canonical choice of parameters in Fig. |l|. We 
compute the spin independent neutralino-proton scattering cross section (used as the 
figure of merit in these experiments), and compare it to projections for the sensitiv- 
ity of Stage 2 detectors (CDMS2||, Edelweiss2[||], CRESST2g, ZEPLIN2|ll) 
which are expected to probe a factor of ~ 5 below the XENON-10 bound. ^ 

We also compare expectations in the HM2DM model with the projected sensitivity 
of the proposed SuperCDMS detector with 25 kg of Ge, and with proposed ton-size 
noble hquid detectors (XENON ||, LUX, WARP||] and CLEAN ||), for which we 
use the sensitivity of Warm Argon Project, with 1400 kg of argon as the benchmark. 

Neutralinos gravitationally trapped in the core of the sun may be indirectly detected 



via their annihilation to neutrinos at neutrino telescopes j56|. Here, we present rates 



for detection of — > /i conversions at Antares|£^ or IceCube|58|. The reference 



experimental level we use is the ultimate sensitivity of IceCube, with a muon energy 
threshold of 50 GeV, corresponding to a flux of about 40 muons per km^ per year. 

3. Indirect detection of neutralinos may also be accomplished by detection of high energy 
gamma rays from neutralino annihilation in the galactic center. Such gamma rays[|59| 
have already been detected by EGRET |^], and will be searched for by the GLAST 
experiment [^]. We evaluate the integrated continuum 7 ray flux above a = 1 GeV 
threshold, and take the GLAST sensitivity of 1.0x10"^'' cm~^s~^ as our benchmark. 

4. Indirect detection of neutralinos is also possible via the detection of anti-particles 
from neutralino annihilations in the galactic halo. Proposed and on-going experi- 
ments include searches for positrons p^] (HEATpS], Pamela] 64] and AMS-02[p3]), 



^For a comparison of the implications of different halo model choices for indirect DM detection rates, 
see e.g. Refs. ^, [li). 

®ln our analysis, we took the 7r-nucleon E term to be 45 MeV. Different values recently suggested by 
other groups can change our predictions by about factor of three |52]. 
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antiprotonsdll (BESSfe^], Pamela, AMS-02) and anti-deuterons (D) (BESSf 
AMS-02, GAPS[|^). For positrons and antiprotons we evaluate the averaged dif- 
ferential antiparticle flux in a projected energy bin centered at a kinetic energy 
of 20 GeV, where we expect an optimal statistics and signal-to-background ratio 
at space-borne antiparticle detectors[45, ItO]. We take the experimental sensitiv- 



ity to be that of the Pamela experiment after three years of data-taking as our 
benchmark. Finally, we evaluate the average differential anti-deuteron flux in the 
0.1 < Tfj < 0.25 GeV range, where Tq stands for the anti-deuteron kinetic energy 
per nucleon, and compare it to the estimated GAPS sensitivity for an ultra-long du- 



ration balloon-borne experiment [g^] (see Ref. |^ for an updated discussion of the 



role of antideuteron searches in DM indirect detection) 



In Fig. 10, we illustrate the various direct and indirect DM detection rates for our 
canonical case with mo = '"11/2 = 300 GeV, Aq = 0, tan/J = 10 and // > 0, where M2 is 
allowed to vary. The M2 value corresponding to the mSUGRA model is denoted by a solid 
black vertical line at r2 = 1, while the HM2DM scenarios for r2 < and r2 > with ^^^h'^ 



within the WMAP range ( |1 . 1\ ) are shown by the green regions on the left and right ends of 
the plot. The MWDM and BWCA solutions are seen as the very narrow green regions near 
r2 ~ ±(0.5 — 0.6). The dotted lines correspond to the sensitivity level for representative 
searches: i.e., the signal is observable only when the model prediction is higher than the 
corresponding dotted line. 

In frame a), we show the spin-independent neutralino-proton scattering cross section. 

We see that for a bino-like neutralino with rriT^ ~ 100 GeV as in the mSUGRA or BWCA 

Z\ 

cases, the signal will only be detectable if the cross section can be probed at the 10~^ pb 
level as envisioned at superCDMS or at 100-1000 kg noble liquid detectors. In contrast, in 
both the IIM2DM regions where the neutralino has a significant higgsino component, the 
cross section is just below the current bound, and should be detectable at CDMS2. This is 
simply a reflection of the well-known result that MHDM has rather large neutralino-proton 
scattering rates, as is typified by the HB/FP region of the mSUGRA model, and further, 
that experimental sensitivity at the 10~^ pb level will probe a wide class of models with 
MHDM (2|]. 

In Fig. we show the expected value of asi{Zip) in the HM2DM model, resulting 
from a scan in niQ and m^/21 keeping tan/? = 10 and Aq = 0, and where we adjust M2 
to get agreement with (|1 . 1| ) , for M2 > with ^ > (upper frame) and M2 < with 
fi < (lower frame). Also shown are the sensitivity limits for CDMS2, superCDMS and 
WARP 1400 kg. The most striking feature of the figure is that the bulk of the points in 
the scan give a cross section around 10"*^ pb, independent of the neutralino mass. This is 
because the increased bino-higgsino mixing necessary to maintain agreement with ( |1.1| ) for 
larger values of compensates for the drop in cross section for larger neutralino masses. 
We also see that the cross sections for negative M2 are systematically lower than those 
for M2 > 0. We have checked that flipping the relative sign between Mi and fj, causes 
a flip in the relative sign between the hZiZi coupling and the HZiZi coupling, so that 
h- and //-mediated amplitudes for neutralino-nucleon scattering interfere constructively 
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nig =300GeV, m^^ =300GeV, tan^=10, Ag =0, \i >0, m, =171.4GeV 
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Figure 10: Rates for direct and indirect detection of neutralino dark matter vs. r2 for toq = 
™i/2 = 300 GeV, with tan/3 = 10, = 0, /i > 0. The various shadings are the same as in Fig. 
The blue curves correspond to the Adiabatically Contracted N03 dark matter halo model, while 
the purple ones are for the Burkert profile. For each experiment, the signal is observable if the rate 
is above the corresponding dotted curve. 



in the positive // case and destructively for negative accounting for the drop in the 
cross section. (Since squarks are heavy, squark-mediated amplitudes are neg ligible.)'^ The 

''In more detail, the flip of the sign of ^ flips the relative sign between the up- and down- higgsino 
components in Zi. Further, since tua ^ A/weak for the HM2DM model, the Higgs mixing angle satisfies 
tana ~ cot so that for tan/3 ~ 10, Reh^ ~ h and Reh'^ ~ //. As a result, there is a flip of the relative 
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cluster of points below m^^ = 100 GeV in the upper frame are points in the white strip at 
mi/2 ~ 0.2 TeV inside the blue region in Fig. ^a). We see that while superCDMS should 
probe all the points for positive values of M2, the increased sensitivity of ton-size noble 
element detectors appears essential in the M2 < case. 

HM2DM: tan^=10, Ag =0, sign(\i)=sign(M2h =171.4 GeV 



SuperCDMS 




WARP 
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Figure 11: The neutralino-proton scattering cross section for the HM2DM model with tan/? 10 
and Aq = 0, but where we scan over mo and TOi/2 for M2 > with fi > (upper frame), and AI2 < 
with /X < (lower frame). For each point on this plot, M2 is adjusted to saturate the observed DM 
density. We retain only those points where fUyj^^ > 103.5 GeV in this plot. 

Turning to indirect dark matter detection, in Fig. ^6), we show the flux of muons 
sign between the ZiZih and ZiZiH couplings. 
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from neutralino pair annihilations in the core of the Sun, again for our canonical choice of 
parameters introduced in Fig. [l|. The blue and purple curves are our projections assuming 
the Adiabatically Contracted N03 Halo model||4^ and the Burkert profile|^], respectively, 
for the distribution of DM in our galaxy. In this case, the result is only mildly sensitive to 
the choice of halo distribution, since the muon flux is mainly determined by the equilibrium 
density of neutralinos in the sun, and the cross section for neutralino annihilation into 
neutrinos. We see that the expected muon flux is more than an order of magnitude below 
the projected sensitivity of IceCube in the mSUGRA framework, but is larger by over two 
orders of magnitude, and in the observable range, for the HM2DM model on account of 
the increased higgsino content of the LSP. The muon flux rapidly drops off near the ends 
where m^^ falls below Mw because the neutrinos (which then mainly come from decays of 
heavy quarks) become very soft, so that their efficiency for detection at IceCube is vastly 
degraded. 

In frames c), d), e) and f) we show the flux of photons, positrons, antiprotons and 
antideuterons, respectively. Also shown by the horizontal lines are the anticipated experi- 
mental sensitivities. Again, we show results using the Adiabatically Contracted N03 Halo 
model (blue) and the Burkert profile (purple). We see that the projections are sensitive 
to the assumed halo distribution. This is most striking for the photon signal at GLAST, 
where the difference is more than four orders of magnitude. This is because most of the 
photons come from the galactic center where the difference between the two profiles is the 
most pronounced. In contrast, projections for the detection of positrons, anti-protons and 
anti-deuterons from neutralino annihilation (unlike photons, these can reach the earth only 
from limited distances) differ by a factor of 5-15. Again, the rates for indirect detection 
via observation of halo annihilation remnants, which are typically low for bino-like DM 
as in the mSUGRA model, jump by a factor of 30-300 in the HM2DM model where \r2\ 
is increased so that the measured CDM relic abundance is obtained. For our choice of 
model parameters, the GAPS experiment should be able to detect anti-deuterons even for 
the case of the Burkert halo profile, while the positron signal in Pamela is projected to be 
just below its sensitivity limit even for the optimistic N03 Halo profile. The situation for 
the anti-proton signal is less conclusive since its detectability clearly depends on the halo 
distribution. 

5. Supersymmetry signals at colliders 

We now turn to an examination of the implications of the HM2DM model for SUSY collider 
searches at the Fermilab Tevatron, the CERN LHC, and a 0.5-1.5 TeV linear e^e^ collider. 
The sparticle mass spectrum in the HM2DM model qualitatively differs from mSUGRA in 
several respects: i) the low parameter implies that charginos and neutralinos should be 
lighter than in mSUGRA cases with gaugino mass unification and a large parameter, 
so these sparticles should be more accessible to collider searches. In addition, ii) the large 
IM2I parameter means W2 and Z4 will be quite heavy and nearly pure wino states, so likely 
difficult to access at colliders, except perhaps via the (kinematically suppressed) production 
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and subsequent decays of ql. Finally, in) the large IM2I parameter pushes left-sfermion 
soft terms to higher values, so that the lighter sfermions are dominantly right-sfermions. 

In Fig. we show contours of in the mo vs. mi/2 plane for ^0 = 0, tan/3 = 10 
and fj. > 0, where at every point M2 has been dialed up to obtain MHDM with ^^_^h'^ ~ 0.1. 
We see that throughout the plane, w-^^ ~ ^^1/27 whereas in mSUGRA, ~ |"ii/2- 

This may be of relevance at an e~^e~ collider where the determination of both and 
^1/2 (via the determination of rn^^) along with M2 may be possible if W1W2 production 
is kinematically accessible. 



HM2DM: > ' tan^=10, =0, \i >0, m, =171.4 GeV 
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Figure 12: Contours of mpj^ in the niQ vs. mi/2 plane for tan/? = 10, Aq = and /i > where 
M2 has been dialed up at every point to yield MHDM where ^^_^h^ ^0.1. 

In addition, in the HM2DM model, since Z2 and Z3 contain large higgsino components 
while Zi is a mixed higgsino-bino state, the mass gaps m^^ — m^^ and — m^^ will 
be expected to be much smaller than in mSUGRA. In mSUGRA, — m^^ ~ 0.4m^/2 
so that as mi/2 grows, the growing Z2 — Z\ mass gap ultimately allows the Z2 — > Z\Z 
or Z\h two-body decays to turn on, which dominate the Z2 branching fraction. The two- 
body "spoiler" decay modes]?^ turn off the leptonic decays Z2 Z\ll^ which can be the 
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starting point for sparticle mass reconstruction in gluino and squark cascade decays [[7^] at 
hadron colliders. In Fig. |l^a), we show the Z2 — Z\ mass gap in the HM2DM model in 
the same plane as in Fig. We see that in the case of HM2DM model, the mass gap is 
everywhere less than M^, so that Z\ll decays will not be shut off by the two-body 

spoiler modes. Moreover, the mass gap is also almost always larger than ~ 25 GeV, and 
decreases with increasing mi/2 in contrast to the cases of MWDM and DM via BWCA, 
where the mass gap increases to beyond 100 GeV for the largest values of mi/2- This could 
serve to distinguish HM2DM from other scenarios, something we will return to below. Also, 
in contrast to mSUGRA as well as the MWDM and BWCA frameworks, in the HM2DM 
model (and other models with MHDM), we expect that m^^ — m^^ cannot be too large 
because the mass of Z3 is expected to be not very far above as we saw in Fig. H In 
Fig. 13&), we show contours of — ra-^^. We see that this difference is also everywhere 
smaller than M^, so that the decays Z3 Z\U and Z3 Z2U are not shut off by the 
spoiler modes. 



HM2DM: > m,,^ , tanf,=10, 4„ =0, \i >0, m, =171.4 GeV 



HM2DM: > m,,^ , tan^=10, A„ =0, n >0, m, =171.4 GeV 




Figure 13: Contours of a) m^^ — m^^, and 6) m^^ — m^^, in the mo vs. mi/2 plane for tan/3 — 10, 
Ao = and > where M2 has been dialed up at every point to yield MHDM where ~ 0.1. 



5.1 Fermilab Tevatron 

In the mSUGRA model, since > 103.5 GeV from LEP2 searches, we expect nig ~ 

3.5m^y^ ~ 350 — 400 GeV, and this high of sparticle masses generally gives quite low 



gg, qq and gq production cross sections|74]. Gluino and right-squark masses are relatively 
unchanged in going from mSUGRA to IIM2DM, while left-squark masses typically increase. 
Thus we expect that gluino and squark production rates at the Tevatron will be relatively 
low, and the signal not easily extracted from data in the HM2DM model. 



Another possibility is to look for pp -^WiZ2^3i + £'f?'""[Q in the IIM2DM model 



The large higgsino components of Wi and Z2 imply that W1Z2 production via the W* 
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will dominantly occur via the isodoublet couplings of the neutralinos to the W. However, 
this (iso-doublet) WW1Z2 coupling will now be smaller than the corresponding coupling 
in the mSUGRA framework (where the coupling arises from the iso-triplet components of 
Wi and Z2, and so is very large) so that a{WiZ2) will be suppressed in HM2DM relative 
to mSUGRA. 

We are thus led to re-examine the rate for trilepton production at the Tevatron in 
HM2DM model and compare this to the mSUGRA case where the trilepton plus EJ^^^^ 
signal is regarded as the gold-plated signature. We use the cuts SC2 proposed in Ref. |75| 
which allow for efficient extraction of the 31 + EJ^^^^ signal from various SM backgrounds, 
the largest of which are W*Z* ^31 + 1^ and W*Y ^ 3£+u production. In Fig. we plot 
the 3i cross section after cuts SC2, fixing mo = 300 GeV, ^0 = 0, tan /3 = 10 and /x > 
versus ^ri^^, obtained by varying mi/2- Contrary to expectation, we see that the signal 
is larger in the HM2DM framework than in the mSUGRA case. We have traced this to 
the fact that the leptonic decay of Z2 is very suppressed within the mSUGRA framework, 
while the corresponding branching fraction is just under B{Z^ ££) ~ 6% as expected 
for the case of MHDM. In the mSUGRA case, the signal rate is always below the 0.8 fb 
level. The total background estimated in Ref. [[f^ j is 1.05 fb, which requires a 1.6 fb signal 
to give a 5a effect with 10 fb~^ of integrated luminosity. The HM2DM curve is enhanced 
relative to mSUGRA, and reaches a maximum of about 1.4 fb - just below the edge of 
observability. The dashed region between rriyy^ ~ 117 — 137 GeV is excluded as can be 
seen from Fig. |^, while for higher rriyy^ values, the signal is always below 0.6 fb, so will be 
difficult to extract at the Tevatron. 



5.2 CERN LHC 

At the CERN LHC, gluino and squark pair production is the dominant sparticle production 



mechanism if sparticle masses are in the TeV range [26|. Gluino and squark cascade decays 
give rise to events containing multiple hard jets, isolated leptons, E^^^^, and sometimes 
also isolated photons. The reach of the LHC has been calculated in the mSUGRA model 



in Refs. |76|. The ultimate reach in the E^^^^ + jets channel is relatively insensitive to the 
details of the cascade decays, but depends mostly on the gluino and squark production 
cross sections, which in turn depends only on their masses. 

We have made a rough translation of the LHC reach in the mSUGRA framework to 
that of the HM2DM model by computing the total production cross section along the 
100 fb~^ reach contour in the last paper of Ref. ||7^, and equating it to total production 
cross sections in the HM2DM tuq vs. mi^2 plane. The result is shown as the dotted contour 



labelled LHC in Fig. 15. At the small mo end, the contour ends up being typically lower 
in mi/2 values by about 150 GeV than the corresponding plots in mSUGRA, with the 
reduction being smaller for larger values of mo. This is because in the HM2DM case, for 
given mo and mi/2 values, the various left-squark masses are raised up by a few hundred 
GeV, causing the corresponding production cross sections to drop. Thus, somewhat lower 
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Figure 14: Trilepton signal after cuts SC2 from Ref. |^ versus TUyy^ along a line of toq = 
300 GeV, variable mi/2, adjusted to reproduce the value of m^^, tan/3 = 10, Aq — and fi > 
where M2 has been dialed up at every point to yield MHDM where ^^i^^ ^ ^-^^ show the 

corresponding rate from mSUGRA model where M2 — mi/2, and the 5cr discovery level for 10 fb^^ 



of integrated luminosity. The dashed region with 117 GeV 
HM2DM model as can be seen from Fig. |^. 



137 GeV is excluded in the 



gluino and right-squark masses are needed to obtain the same reach level as mSUGRA 
within the HM2DM framework. The location of the LHC reach contom' in the HM2DM 
model varies between mi/2 ~ 1260 GeV at low mo and ~ 1040 GeV at mg = 2 TeV. 
This corresponds to a reach in terms of gluino mass of rrig ~ 2350 — 2750 GeV for 100 fb~^ 
of integrated luminosity. 

If supersymmetry is discovered at the LHC, reconstruction of SUSY events to measure 
sparticle masses will be an important item on the experimental agenda. Edges in the 
distribution of opposite sign, same flavour dilepton masses provide important information 



on the mass difference m? 



m-r 



between those neutralino pairs for which the three body 



decay Zi ZjU has a significant branching fraction |77]. In the mSUGRA framework, 
we typically expect an observable edge just for the decays of Z2, if its mass is low enough 
so that the spoiler two body modes are not accessible. However, for the HM2DM model, 
as for all models with MHDM, both Z2 and Z^ are relatively light, so that we may expect 
more than one mass edge in the dilepton mass spectrum. 

To examine this, we simulated IM LHC SUSY events for the HM2DM1 point in Ta- 
ble H using Isajet 7.76, and passed these through a toy detector simulation as described 
in Ref. |78|. In addition to the various geometric and other acceptance requirements, we 
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HM2DM: > m^^ ' tan^=10, Ag =0, ^ >0, m, =171.4 GeV 
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Figure 15: Approximate reach of the CERN LHC (with 100 fb"^ of data) and ILC with ^/s = 0.5 
and 1 TeV in the HM2DM model, viewed in the mo vs. toi/2 plane for tan/3 = 10, Aq = and 
fi > where M2 has been dialed up at every point to yield MHDM where ^^i^^ ~ 0.1. 



required the following analysis cuts on the signal: 

E:^''' > (100 GeV,0.2Meff), 

"-jets > 4, 

ETiji,j2,j3,k) > (150,150,80,50) GeV, 
St > 0.2. 

For the present analysis, we only retain events which include, in addition, a pair of opposite 
sign (OS) dileptons. The invariant mass distribution for e~^e~ + /^"^/U" pairs in SUSY 
events is shown by the red-shaded histogram for the signal, and by the line histogram 



for the SM background, in Fig. 16 a). In the signal, these leptons arise from the decays 
of neutralinos and charginos produced in the SUSY cascade. The mass of dileptons from 
the three body decay Zi Zj££ is kinematically bounded by — , resulting in the 
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mass edges mentioned above. Lepton pairs where the leptons each come from decays of 
different charginos or neutrahnos are uncorrelated in mass, and so are expected to yield a 
smooth broad continuum. Since these pairs are also uncorrelated in flavour, subtracting the 
distribution of unlike- flavour OS dilepton pairs should (up to statistical fluctuations) leave 
us with the dilepton mass distribution of dileptons from the same neutralino ||7^]. This 



subtracted distribution is shown in Fig. 166). The following features are worth noting. 



1. The red shaded histogram in frame a) shows distinct gaps near around 46 GeV and 
at around 64 GeV, close to the mass edges at 47.6 GeV and 65.4 GeV expected from 
the leptonic decays of Z2 and Z3 to the LSP, respectively. There is no corresponding 
gap near m^^ — m^^ indicating that the branching fraction B{Z3 Z2U) is very 
small (0.2% per lepton flavour in the present case). 

2. There is a distinct peak at Mz showing that Z bosons are being produced in SUSY 
cascades. This is a strong indication (even if not evidence) of the production of W2 
and/or Z4, since the only other source of Z's would be the decays ^2 tiZ, b2 ^ biZ 
or ^4 — > hZ, all of which are likely to have very small cross sections in many models. 

3. The subtracted distribution in frame b) also shows these edges even though the 
location of the edge seems somewhat smeared out by an upward fluctuation 
in the opposite flavour OS dilepton distribution in the mass bin just below 60 GeV, 
resulting instead in a shoulder around the expected value. 

4. The large levels of the background in some bins are caused by a handful of QCD 
events passing our cuts. Because of the large QCD cross section, we need to simulate 
a much larger sample of QCD events than we were able to in order to further reduce 
these fluctuations.^ We expect though that the QCD background to the SUSY signal 
will be under control at te LHC. 

5. Finally, we observe that in both frames the shape of the dilepton distribution in the 
region below the first edge, and in between the two edges is very different. Specifically, 
the — m^^ mass edge is much sharper than the m^^ — edge. This is due to 
the difference in shapes of the mass distributions of dileptons from the decays of Z2 
and Z3. As pointed out by Kitano and Nomura in Ref. [^0|, this shape depends on 
the relative sign of the mass eigenvalue of the parent and daughter neutralino. If this 
relative sign is positive, we get a sharp mass edge as is the case for the ^z3~^Zi ^"^S^ 
in the figure. Since Z2 and Z3 are dominantly the higgsino states, not surprisingly 
they then have opposite signs for the eigenvalues, causing the distribution of dileptons 
from Z2 iiZi to peak at lower values of m(£~^i~) so that the m^^ ~ ''^Zi ^^S^ 
much more diffuse. The shape of the dilepton distribution thus provides important 
information about the composition of neutralinos. 



*It is not clear that such a computer-intensive simulation of the tail of this detector-dependent back- 
ground would be particularly meaningful, especially in view of our toy calorimeter simulation. 



5.3 Linear e^e collider 



In assessing the role of the ILC for a discovery of SUSY within the HM2DM framework, we 
first note that, for a fixed value of mo and the value of fn^^ is lowered with respect to 

mSUGRA. This is because in the HM2DM model, is lowered to below ~ 2Mi, which is 
approximately the Wi mass in models with unified gaugino masses. This increases the reach 
of ILC in mo vs. mi/2 space while at the same time the LHC reach is slightly diminished 



due to the increased left-squark masses. We show in Fig. 15 the approximate ILC reach 
for a y/s = 0.5 and 1 TeV machines by delineating the mass contours where and 
nif-^ = 250 and 500 GeV (some additional region may be accessible beyond this via Z1Z2 
production |8^). The bulk of the ILC reach for a ^/s = 0.5 TeV machine reaches to the 
1^1/2 ~ 500 GeV level, corresponding to a gluino mass of ~ 1150 GeV. A ^/s = 1 TeV ILC 
has a reach extending to mi/2 ~ 1150 GeV, corresponding to a value of nig ~ 2600 GeV. 
The enhanced reach of ILC coupled to decreased reach of LHC in terms of mi/2 means 
that a 1 TeV ILC has a comparable reach to the LHC. 

Within the HM2DM model, we have Mi ~ < M2 at the weak scale. Thus if 
gluinos ( remember that rrig ~ 6M1) are accessible at the LHC, it is reasonable to expect 
that Wi as well as Zi, Z2 and Z3 will be accessible at a TeV linear collider. The HM2DM 
framework will be readily distinguishable from the mSUGRA model if chargino pair pro- 
duction is kinematically accessible. Since the chargino is mainly a wino in mSUGRA, and 
a higgsino in the HM2DM model, the total chargino pair production cross section (for a 
given value of 'm^^^ as well as its dependence on the polarization of the electron beam is 
very different. In the mSUGRA case, the expected cross section is significantly larger than 
in the HM2DM model, and further, drops to very low values as the electron beam is taken 
to be increasingly right-handed, while this same dependence is comparably milder for the 
hieesino-like W^ . This mav be corroborated by studying the polarization-dependence of 




2, ^ ^1 





m„ (GeV) 



Figure 16: a) Distribution of the invariant mass of same flavor /opposite sign dilepton pairs in 
SUSY events at the CERN LHC (red histogram) after the cuts discussed in the text, along with 
the expectation for the corresponding SM background (open histogram) . b) The difference between 
this distribution and the corresponding one with opposite flavour/opposite sign deleptons in SUSY 
events (blue histogram) along with the SM background (open histogram) . The signal distributions 
are calculated for the HM2DM1 case in Table 0. 
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neutralino pair production. Within mSUGRA, the polarization dependence of the wino- 
Uke Z2 pairs is similar to that of chargino pair production, while the higgsino-like Z3 and 
Zi^ are typically heavy. For the MHDM case realized in the HM2DM framework, on the 
other hand, Zi, Z2 and Z3 are light and mixed, so that not only are several ZiZj pairs 
accessible, the polarization dependence and size of the cross sections is very different. 



This is illustrated in Fig. 17, where we show various -ino pair production cross sections 
accessible to a -y/i = 0.5 TeV ILC versus beam polarization PL{e~) in the case of a) the 
mSUGRA model case and b) for the HM2DM1 case in Table |I[ In frame a), we see that 
only W^W^ , Z1Z2 and Z2Z2 are accessible to a 0.5 TeV ILC, and that their cross sections 
precipitously drop from readily observable values of tens or hundreds of femtobarns to below 
1 fb, as the electron beam polarization becomes increasingly right-handed. In the IIM2DM 
case shown in frame b) , we see that as anticipated many more -ino pair production reactions 
are accessible, with vastly differing cross sections. The production of one chargino and 
three neutralino states should be unambiguous. The relative size of the various neutralino 
cross sections will be parameter-dependent: in the present case the small size of (7(^2^2) 
is a reflection of the fact that the magnitudes of the and components of Z2 are 
nearly equal. The polarization dependence of the (higgsino-like) neutralino production 
cross sections is also different from frame a). A detailed study of the -ino production 
reactions should allow the determination of Mi and 



mSUGRA 

=300, m|,,=300, Aj^^O. tanP^lO, .sgnn=+ 



HM2DM 

=300. m =300, A =0, IimP=!0, sgnfi=+l; M,=900 





Figure 17: Cross section for W^W^ and ZiZj production at a \fs = 0.5 TeV linear collider versus 
beam polarization parameter P^le ) for toq — 300 GeV, mi/2 = 300 GeV, — 0, tan/3 — 10 and 
/,t > in the a) mSUGRA model and 6) HM2DM model with M2 = 900 GeV. 

Another feature of the HM2DM model is that, since M2 is large and feeds into sfermion 
masses via the RG running from a universal scalar mass, we expect the lightest sfermions 
to be dominantly right-type, and the heaviest are dominantly left-type. Within mSUGRA, 
the sfermions of the first two generations are approximately degenerate, while for the 
third generation only the r-sleptons and top squarks are dominantly right-handed, while 
large top quark Yukawa effects make b\ mostly the left-type, as illustrated in Fig. |6|. We 
illustrate the beam polarization dependence of third generation sfermion pair production 
cross sections at a hypothetical -y/i = 1.5 TeV ILC in Fig. |l^ for a) the mSUGRA model 
point in Table || and 6) the case HM2DM2 case in Table |. It is clear in the mSUGRA 
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model from the polarization dependence that in fact fi is dominantly right-type, bi is 
dominantly left-type and ti is mixed left-right. While the polarization dependence for stau 
and stop pair production is qualitatively similar in the two models, that for 6-squark pair 
production is markedly different, providing a clear signature for the qualitatively different 
value of 9b in the HM2DM model. 



mSUGRA HM2DM 

m^-300, nij^^-SOO, Aj^-O, lanP-10. sgn^-+l m^-300. m|^,-300. Ajj-O. tiinP-10. sgn(l-+l; M,-900 




P P 

L L 

Figure 18: Cross section for iiii, and fif{' production at a ^/s = 1500 GeV linear collider 
versus beam polarization parameter PL{e~) for mo = 300 GeV, mi/2 — 300 GeV, Ao — 0, tan/3 — 
10 and > in the a) mSUGRA model and b) HM2DM model case with M2 = 900 GeV in Table 0. 



6. Summary and concluding remarks 

While the recently measured value of the CDM relic density can be accommodated within 
the paradigm mSUGRA framework, much of the allowed region lies at the edge of the 
parameter space as exemplified by many studies in the mo — ^1-1/2 plane. There are, 
however, several one-parameter extensions of the mSUGRA model where it is possible to 
get the observed value of Ocdm^^ all over this plane for essentially all values of tan (3 and Aq. 
These extensions involve either the adjustment of the mass spectrum (so that the neutralino 
annihilation rate is enhanced via an s-channel A/H resonance or via co-annihilation with 
a charged sparticle) or the adjustment of the neutralino composition (to obtain either 
MWDM or MHDM). MHDM requires a reduced value of and several models that lead 
to small 1^1 have been proposed. In this paper, we have pointed out a new mechanism for 
obtaining the observed relic density: non-universal boundary conditions with a large GUT 
scale value of the SU{2) gaugino mass parameter IM2I can lead to MHDM. We have also 
studied the broad phenomenological implications of this scenario. 

Common to all scenarios with MHDM, our scenario has an -ino spectrum where the 
lightest neutralino is a mixed bino-higgsino state, the lighter chargino and the next two 
heavier neutralinos have large higgsino content, while the heaviest chargino and neutralino 
are dominantly wino-like with a mass IM2I considerably larger than the other -ino masses. 
In the HM2DM model, IM2I is raised even further and the wino states become very heavy 
and would likely only be produced at the LHC via cascade decays of ql- The feature that 
distinguishes the HM2DM model from other models with MHDM is that weak interaction 
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effects make the left-sfermions significantly heavier than their right-sibhngs. This effect 
leads to a qualitative change in the intrageneration mixing pattern of 6-squarks: in the 
HM2DM model, 61 is dominantly right-handed, while in most models (with universal mass 
parameters for and ft^j) top-quark Yukawa coupling effects make 61 mostly left-handed. 

As in all models with MHDM, direct search experiments provide a promising av- 
enue once they reach a sensitivity to probe neutralino-nucleon scattering cross sections 
~ 10~® pb. With an order of magnitude increase in sensitivity (projected at superCDMS 
or at ~ 100 kg noble liquid detectors), there should be an observable signal over most 
of the parameter space, while ton-sized noble element detectors should be able to probe 
the entire parameter space. Our conclusions for indirect searches are less definitive since, 
except for IceCube type detectors, the signals depend on assumptions of the distribution 
of the DM in our galactic halo. While IceCube may well be sensitive to the signal, the best 
prospects for anti-particle detection appear to come from the search for anti-deuterons. 
For the case that we examined, the Pamela satellite that is currently gathering data may 
just be sensitive enough to the signal from anti-protons for a favourable halo profile, but 
less so for the positron signal. The dependence of the gamma ray signal from our galactic 
center on the halo profile is too large to draw strong conclusions about the observability 
in GLAST, but perhaps a signal in anti-particle searches may make the situation clearer. 

We have discussed collider signals in Sec. |5|. Prospect for a discovery at the Tevatron 
are not encouraging within this scenario. While the reach of the LHC is slightly degraded 
relative to its reach in the mSUGRA model, experiments at the LHC should be sensitive 
to gluino masses as large as 2750 GeV (2350 GeV) for a small (moderate) value of mo, 
assuming an integrated luminosity of 100 fb. The TeV linear collider should be sensitive to 
a discovery in a parameter space region very similar to the LHC. Searches in the multi-jet 
plus opposite sign dileptons plus EJ^^^^^ channel will be especially interesting as the dilepton 
mass distribution for same fiavour lepton pairs may not only allow the construction of more 
than one mass edge (strongly suggestive of MHDM), but may also, via its shape, provide 
indirect evidence for a higgsino-like Z2/Z3. 

It should be possible, at least in principle, to distinguish the HM2DM model discussed 
here from other scenarios that also lead to agreement with the observed relic density. 
Within a supersymmetric interpretation of an observed signal at the LHC, the observation 
of more than one dilepton mass edge would point to the existence of Z2 and Z3 with 
a relatively small mass splitting (< Mz) between Zi and these neutralinos. In turn, this 
would point to a small value of If we further assume that the CDM density is saturated 
by the LSP, we know that Zi cannot be the higgsino. Recall that for Mi <C <^ M2, 
m^^ ~ Ml while and would be on either side of \n\ so that the value of — 
must lie between the — m^^ and — mass edges. We note, however, that 
while the second inequality always holds within the HM2DM model, |//| may not always 
be that much larger than Mi as, for instance, in the HM2DM1 case. The mass edges 
nevertheless serve as a semi-quantitative indicator of |;u| relative tom^^, and so may provide 

^It is, of course, logically possible that M2 is just a bit smaller than so that all four neutralinos are 
strongly mixed and accessible, but with m^^ — m^_^ > Mz- In this case it would not be unreasonable to 
expect additional mass edges from decays of the heavier neutralinos to Z2, or perhaps even to Z3. 
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corroborative evidence for consistency with the MHDM hypothesis. A striking confirmation 
of the MHDM hypothesis could come from the observation of a large cross section in direct 
DM search experiments. While the dilepton mass edges may serve to separate out the 
MHDM scenarios from those with MWDM or BWCA, they do not provide evidence for 
the HM2DM framework, since this would likely require a complete reconstruction of the 
-ino sector with a precision difficult at the LHC. Aside from a determination of M2 via that 
of (or fn^^, the smoking gun for the HM2DM framework would be observing a wino- 
like W2/Z4 and/or a large splitting between the left and right sfermions, or determining 
that the lighter 6-squark state is mainly bR. None of these appear straightforward at the 
LHC. It is unlikely the heavier slepton or chargino will be accessible at the LHC in this 
scenario, and squark mass splittings and mixing angles are difficult to measure there. A 
high energy electron-positron collider offers the best prospects for these measurements, 
assuming, of course, that these sparticles are kinematically accessible. 
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